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these problems, : 

This thesis investigates a model of a slot antenna using the Numerical Elec- 
bronagnetics Code(NEC). Input impedance, near field inside the slot and ra- 
ciation patterns are presented, 
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ABSTRACT 


One of the major problems in communications aboard ships 
is the minimization of the number of antennas, avoiding 
undesirable coupling between them and interference between 
links, therefore reducing the exposed area of ships. 

In the future, for the next generation of ships the 
trend will be the elimination of tall and large structures 
to make antennas more survivable during combat. 

The use of slot antennas along the bulkheads of the 
superstructure апа broadband techniques applied to them 
seems to be a good candidate for solving these problens. 

This thesis investigates a model of a slot antenna using 
the Numerical Electromagnetics Code (NEC). Input impedance, 
near field inside the slot and radiation patterns are 


presented . 
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I. INTRODUCTION 


A. SCOPE OF THE STUDY 

One of the major problems in communications aboard ships 
is the minimization of the number of antennas, avoiding 
undesirable coupling between them and interference between 
links. 

In the areas of antenna and propagation, however, major 
advances in performance have not been easily achieved. This 
is due mainly to the non-deterministic behavior of 
ionosphere and because of antennas size limitations imposed 
by survivability and/or visibility profile requirements. Le 
should not be forgotten that electrically small antennas 
have narrow bandwidth, strong platform interaction and 
undesirable coupling to cosited antennas. 

In the last years, the decks and superstructures aboard 
naval ships have become much too cluttered with antennas of 
all types and sizes. It can be said that, in general, one 
antenna was installed for each HF transmitter, transceiver, 
or receiver. 

Major advances in communication equipment performance 
and size reduction have been achieved because of solid state 
technology. But efficient antennas have to be on the order 
of a quarter wavelength in size and antenna miniaturization 
is not easily achieved at the longer wavelengths. (At a 
frequency of 2 MHz, a quarter wavelength vertical antenna 
would have to be 37.5 meters high. Such antenna might be 
seen as "elctrically small" but not small physically). 

Mutual coupling of cosited antennas can deteriorate the 


performance of communication systems involving a number of 


radios operating on adjacent frequencies. The situation 
ο οἴ οἱ δα in shipborne platforms where simultaneous 
transmission and reception takes place. Crosstalk 
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interference due to mutual coupling of cosited antennas is 
often very difficult to mitigate in such installations. 

The elimination of tall and large structures will make 
antennas more survivable during combat. 

Im the future, for the next generation of ships, the 
trend will be to build ships without masts, using their 
bulkheads along the superstructure to deploy, for example, 
slot antennas which seems to Бе good candidates for solving 
these problems. 

This paper is an investigation of the performance of an 
HE slot deployed on a simulated bulkhead, created by a 
Wire-grid wall, using a special computer program called 
NEC-Numerical Electromagnetics Code. (A brief description 
of NEC is given in Chapter III). 

The understanding of linear antennas has been advanced 
since the late 1960's by the introduction of the METHOD OF 
MOMENTS. Versatile computer codes based on this method, 
like NEC, allow very complex antenna systems to be analyzed. 

In Chapter V results for input impedance, radiation 
pattern, and bandwidth of the slot antennas are shown. Ап 
attempt to find the best way to feed the slots using the 


standard transmission lines is another goal of this thesis. 


B. REVIEW OF PAST WORKS 

In 1954, A. Young studied [Ref. 1] how to solve integral 
equations using approximate product-integration, without the 
aid of a digital computer. 

K. K. Mei published a paper [Ref. 2] in 1965 describing 
the use of Hallen's integral equation for thin wire 
antennas, during a time when large digital computers were 
emerging. 

To solve scattering problems using Pocklington's 
integral equation , one procedure is to reduce the integral 
equation to a system of simultaneous linear equations. R. F. 
Harrington in 1967 used the Method of Moments [Ref. 3] for 


JE IE 


treating field problems. He described the use of pulse and 
piecewise linear current basis functions to solve the scalar 
and vector potential functions. 

In 1950, М. Begovich studied the resonant slot antenna 
in an infinite plane conductor [Ref. 4]. 

In 1978, G.J.Burke and others, prepared a paper [Ref. 5] 
describing a computer program entitled Numerical 
Electromagnetics Code for analyZing the electromagnetic 
response of an arbitrary structure consisting of wires and 
surfaces in free space or over a ground plane. 

Before the advent of the computer, electromagnetic 
problems were formulated in terms of partial differential 
equations with boundary conditions or in terms of an 
integral equation. Much work was needed to apply these 
classical methods to practical problems and only antennas 
with simple geometries could be treated. 

With the Method of Moments, the integral equation for 
the antenna current is replaced with a matrix equation and 
then solved. Computer codes such as NEC and MININEC now 
enable the analysis of very complicated antenna systems and 
provide answers to a good degree of accuracy. NEC requires a 
main-frame computer for efficient implementation. MININEC 
was written in BASIC for PC's. Тһе МЕС апа МІМІМЕС have 
taken the electromagnetic community "by storm" because of 
their great versatility and applicability. With these codes 
the antenna enginneer can test different approaches to 
design problems and can explore new trends and concentrate 
on his immediate problems without getting involved with the 
complications of classical elctromagnetic theory. A good 
fundamental knowledge of electromagnetics, some familiarity 
with methods and limitations of the antenna computer code, 
and a certain experience with NEC for known cases, including 
the ability to validate answers from the code, are important 


requirements for the engineer. 


T2 


Numerical Electromagnetics Code (NEC)-Method of Moments 
is a paper of G.J.Burke and A.J.Pogio [Ref. 7] published in 
1981, describing a user-oriented computer code for analysis 
of the electromagnetic response of antennas and other metal 


surfaces. 
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Il. THEORETICAL BACKGROUND 


А. SLOT ANTENNAS 
1. Zhe Principle of Slot Antennas 

Slot antennas may be cut in a metallic surface such 
as the bulkhead of a ship, the skin of an aircraft or the 
wall of a waveguide. The slot may be fed by a transmission 
line or a generator connected across it. In this paper a 
single generator is used to feed the slot antenna. 

The antenna built with two resonant stubs of quarter 
wavelength connected to a 2-wire transmission line (Fig.2.1) 
forms a poor radiator. If the lines are spaced very close 
and carry currents of opposite phase, their fields cancel 
one each other. The end wires carry currents in same phase 


put they are too short to radiate efficiently. 





Figure 2.1 Poor Radiator. 


But now the antenna shown in Fig.2.2 сап Бе 


considered as a very efficient radiator. 


14 


Metal sheet 


y 





Figure 2.2 Efficient Slot Radiator. 


This is a slot antenna. The currents are spread out 
all over the conducting surface and there is radiation from 
both sides. For horizontal slots the radiation 
perpendicular to the sheet is vertically polarized. The 
input impedance is about 500 ohms for a resonant slot of 
half wavelength fed at its center. It 1s possible to 
decrease the input impedance with an off-center feed point. 
For a coaxial cable of 50 ohms this offset distance should 
be about S-1/20 of the wavelength as seen in Fig. 2.3. 

As was stated before, the conducting sheet radiates 
equally from both sides, but if a very large boxed sheet is 
placed at one side of the slot, the radiation will occur 
only from one side (see Fig.2.4). 

If the depth of the box placed behind the slot is 
about one fourth wavelength, there will be no susceptance at 
the terminal of the slot and the slot input impedance should 
be 1000 ohms [Ref. 8]. 
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Figure 2.3 Special Feeding. 





is 
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the 


Figure 2.4 Boxed Slot. 


Radiation Pattern of Slots 
Booker [Ref. 9] showed that a slot radiation pattern 


same as that of а complementary horizontal 
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half-wavelength dipole (see Fig. 2 5) Consisting of а 
perfectly conducting flat strip of same width. 
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Figure 2.5 Complementary Pair. 


But there are two differences: 

e the electric and magnetic fields are interchanged and 
EMEN εξ Έτοποπα side to the 
EM αρα magnetic field. — e 

Ihe infinite flat sheet is coincident with x-z plane 
and the long dimension of the slot is in the x-direction 
(see Fig. 2.6). 

The radiation-field pattern has the same doughnut 
shape but the direction of E, for example, is interchanged 
if compared with that of the complementary dipole antenna 
coincident with x-axis. 

J. Ihe Impedance of Slot Antennas 

By using Babinet's principle, many of the problems 
of slot antennas can be reduced to situations involving 
complementary linear antennas for which solutions have 
already been obtained. 


% 


yy 
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Figure 2.6 E-Field Radiation Pattern. 


Kraus showed how to develop a relation for the 
impedance Zs of a slot antenna in terms of the impedance Zd 
of a complementary dipole antenna. 


Since Zd may be complex: 


Zs = 20° /42d = 376.72/4(Rd + jXd ) = (еап 2210 


35,476(Rd - jXd )/( Rd? + Xd? ) 


It can be easily seen that if the dipole antenna is 
inductive, the slot is capacitive and vice-versa. 

Lengthening a half wavelength dipole makes it more 
inductive, but lengthening a half wavelength slot makes it 
more capacitive. 

Kraus gives some numerical examples that will Бе 
used later on comparing with NEC results: 


e the terminal impedance of an infinitesimally thin half 
wavelength slot antenna 45: 


Zs =35,476/(73 + j42.5)=363 - j211 ohms (eqn 2.2) 


e the terminal resistance of che complementary slot 
antenna compared to that of a cylindrical antenna with 
length-to-diameter (L/D). ratio of 100 and length of 
about 0.475 wavelengths is: 


Zs =35,476/67 = 530 ohms (eqn 2.3) 


uice terminal resistance of a complementary slot antenna 
compared to that of a cylindrical dipole with L/D=28 
and length of about 0.925 wavelengths (Zs = 710 ohms 
is about 50 ohms, which could be easily matched to a 5 
ohms-coaxial cable. 


Finally it must be said that: 


LU ο ο. ους enclosed on one side of the sheet with 
a box of such size that zero  suceptance is shunted 
across the slot terminals, due to the box, the 
impedances are doubled; 


e the bandwidth (BW) or selectivity of a slot antenna are 
the same as for a complementary dipole. If a slot is 
widened (smaller L/W ratio), the slot antenna BW is 
increased the same as increasing the thickness of а 
ng antenna (smaller L/D ratio) increases its BW; 
and, 

e this study applies to slots in infinite sheets. The 


impedance values are the same provided that the edge of 
the sheet is at least a wavelength from the slot. 


B. PRELIMINARY CONSIDERATIONS: BOUNDARY CONDITIONS 
It is worthwhile to review some ideas about boundary 
conditions which are considered by NEC when dealing with 


wire structures. 
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Two main considerations are assumed in this study: 


° the flat surface of a real bulkhead was modeled by the 
wire grid; and 


e the slot opening was outlined by wires. 

These assumptions produce, for example, near fields 
inside the slot and currents on the wire grid which can be 
an acceptable approximation to the real slot placed ina 
solid planar surface. The differences are in the nature of 
concentration of currents along the wires of the grid as 
opposed to a smoother surface current distribution of a 
plate. The near field plot shown in Chapter V is an evidence 
or that Басе. 


20 


IRA BRIER DESCRIPTION OF NEC 


Numerical Electromagnetics Code, NEC, is a computer 
program that analyzes the electromagnetic response of 
antennas. It computes a numerical solution to integral 
equations that describe the currents induced on a structure 
by voltage or current generators or incident fields. 

NEC was developed by Lawrence Livermore Laboratories 


with funding from the Naval Electronics Systems Command. 
A. MAIN CHARACTERISTICS OF NEC 


An Integral Equation (I.E.) for current modeling on thin 
wires is combined with an I.E. for modeling the current on 
smooth surfaces to describe the electromagnetic response of 
an arbitrary structure. The structure may have active parts, 
may be located over a perfect or imperfect ground plane, and 
may have distributed or lumped-element loading. Тһе 
excitation may come from voltage or current sources on the 
structure, an incident plane wave with linear or elliptic 
polarization, or the field may be due to a elementary 
dipole. The outputs may include current and charge density, 
power gain or directive gain, near or far electric or 
magnetic fields, impedance or admittance, total radiated 
power or input power. 

NEC utilizes the Gauss-Doolittle Method for solving the 
matrix equations generated by the Method of Moments when 
solving the I.E.'s. 

It is a discrete sampling code where a complex structure 
must be dissected into a number of simple elements, which 
are thin wires or small surface areas, to which the 
Electric Field Integral Equation (EFIE) or the Magnetic 
Field Integral Equation (MFIE) are applied, respectively 


when using the approximations. The resemblance with the real 
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world is strongly influenced Бу the choice of the 
discretization of the structure by the user. Smaller 
geometric elements produce models which come close to 
reality. But, the smaller the elements, the larger the 
number of elements will be, which means the larger the 
matrix of equations and hence the more expensive the 
solution will be. 

The choice of proper dissecting is gained with 
experience and becomes an art as well as a science. 


Generally, segment lengths (A) should be less than 0.1 


wavelength [Ref. 6]. Shorter segments (0.05 wavelength ) 
may be needed at curves or junctions.  Segments smaller than 
O. 001 wavelength should be avoided. Ihe radius of the wire 


"a" relative to the wavelength depends on the Kernel used in 


the I.E.. There are two options. The thin wire kernel models 
а filament current, while the extended Kernel models а 
uniform current distribution around the segment surface. The 
field of the distributed current is approximated by the 
first two terms in a series expansion of the exact field, in 
powers of a“. The first term (a?) is identical to the thin 
wire kernel; the second term extends the accuracy for larger 
values of a. Both kernels incorporate the thin wire 
approximations and both also require that the ratio between 
the wire perimeter and the wavelength be much less than 1. 
The thin wire kernel requires a ratio /Ma»2. lhese values 
ensure errors are less than 17. It must be said that the 
extended kernel also permits a value of A /a>2. 

The conducting surfaces [Ref. 10] are modeled by small 
flat surface patches which conform as closely as possible to 
curved surfaces. The parameter defining a patch is a normal 
unit vector, originating from the center of the patch, 
defined in cartesian coordinates. Each patch which connects 
to a wire have the connection point at the patch center. 


The code divides each patch into 4 equal patches about the 
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wire end, along the unit vector lines describing the surface 
of the patch. An interpolation function is applied to the 4 
patches to represent placement of the current onto the 
patches and the function is numerically integrated. Patches 
with wires connected to them should be chosen to be 
approximately square with sides parallel to the unit vectors 
defining the surface. Only one wire may connect to a patch 
and it may not lie in the plane of the patch. A minimum of 
about 25 patches should be used per wavelength of surface 
area; the maximum size for an individual patch is about 0.04 
square wavelengths. Long narrow patches should be avoided. 

For a perfectly conducting ground, the code generates a 
reflected image. Structures may be close to, or contact, the 
ground plane. 

The EFIE is used for thin wire structures with small 
conductor volume. It is derived from the electric field 
representation for a current distribution confined to the 
surface of a perfectly conducting body and is solved in NEC. 
The MFIE is derived from the integral representation for 
the magnetic field of a closed surface current distribution 


and is also solved in NEC. 


B. METHOD OF MOMENTS 

Ihe Method of Moments is a technique whereby an integral 
or differential equation is reduced to a system of linear 
algebraic equations which are easily manipulated by high 
speed digital computers. Therefore its objective is the 
conversion of a continuous problem into a discrete problem 
having a finite number of unknowns. 

References [10] and [11] discuss the method showing that 
the accuracy of solution depends on the choice of basis 
function or, in "other words, upon both resolution of the 
discretizement (no. of unknowns) and also upon the method of 
approximating the continuous solution and how the equation 


is enforced. To generate the system of linear equations, 
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the structure is divided into a number of segments. The 
current distribution on the antenna is approximated as the 
summation of currents found on each segment. The next step 
is the identification of a set of weighting functions. These 
are set equal to the basis functions, when using Galerkin's 
Method. In NEC weighting functions are chosen as a set of 
Dirac Delta functions for efficiency of solution. The 
system of equations is then solved for the exact value at 
match points anda very close approximation to the actual 
current distribution at nearby points can be obtained. 

After determining the currents the other non-zero scalar 


field components can be computed. 


C. THIN WIRE APPROXIMATION 

For a thin perfect electrical conductor cylinder with 
radius "a" and length "L", where a << one wavelength and a 
<< L, апа in presence of a known incident field it is 
possible to assume that induced vector surface current, Js , 
is only z-directed and resides on the surface. 

The following thin wire approximations were used: 


° transverse currents can be neglected relative to axial 
currents on the wire; 


• the current can be approximated by filaments at the 
geometrical center the wire axis; 


е electrical field Conese ЕСІН both inc ide ma and 
scattered fields. As the induction theorem states that 
the induced surface current density must producer 
electric field whose tangential component cancels that 
of the incident electric field, on the surface of th 
. the boundary condition on the electric tiem 
to be inforced in the axial direction only;and 

e the circumferential variation in the axial current can 


be neglected and the current at the епаз о1 суллавштин 
surfaces 1s approximately zero. 


D. NUMERICAL GREEN'S FUNCTION 

The Numerical Green's Function (NGF) option is used in 
this paper. A fixed structure and its environment are 
modeled and the factored interaction matrix is saved ina 


file. New parts are added to the model in subsequent 
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computer runs and the complete solution obtained without 
repeating calculations for the data on file. The main 
purpose of the NGF is to avoid unnecessary repetition of 
calculations when a part of the model, such as the 
simulation of the bulkhead and the slot using the wire grid 
model, is modified just by adding a feed segment one or more 
times while the environment remains constant. The 
self-interation matrix for the fixed environment may be 
computed, factored for solution, and saved on a tape or disk 
file. Solution for a new feed point, for example, then 
requires only the evaluation of the self-interation matrix 
for the fed segment, the mutual antenna-to-bulkhead 
interations, and matrix manipulations for a 
partitioned-matrix solution. 

Another reason for using the NGF option is to exploit 
partial symmetry in a structure, which was done in this 
paper. Such partial symmetry may be exploited to reduce 
solution time by running the symmetric part of the model 
first and writing a NGF file. The unsymmetric parts may then 
be added in a second run. 

The use of NGF in this paper during the simulations of 
large, time-consuming models saved expensive results for 
further use. Without adding new wires they were used with 
new excitation, ок to compute new radiation patterns, or 


near field. 
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IV. MODEL DESCRIPIION 


Using NEC geometry cards, the first step of this 
investigation was to build, with wires, a bulkhead of a ship 


where a slot antenna was placed. 




























THETA = 60.00 PHI= 60.00 ETA = 90.00 


Figure 4.1 Bulkhead with Slot Antenna. 
“ 
Thus, the wire-grid modeling ability of NEC is used to 


simulate the behavior of a rectangular slot antenna. 
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ETE. 4.1 shows the bulkhead wire-grid model with 
dimensions 12m x 9.2m in the center of which was placed a 
slot with dimensions 7m x 0. Zm. These are typical bulkhead 
dimensions for river patrol boats or small corvettes. A 
simple rectangular slot cut into a square conducting sheet 
was chosen as the antenna structure to be evaluated using 
the wire-grid modeling approach. This geometry simplifies 
the modeling procedure and , as shown in Chapter II , there 
is some qualitative information available regarding the 


characteristics of rectangular slot antennas. 


BULKHEAD WITH SLOT ANTENNA 

1.1» -6.0,6.6» -%.0.3.6» .01 
2.1» -6.0.6.%. -5.0.4.%6. .01 
2,11, 0.0.0. 1.0.0 

25.1. 6.0.6.6, 6.0.3.6, .01 
25.2. 0.0.0. 0.0,-1 

76.1. -%6.0,1.6» -6.0.,1.1. .01 
77.1» -%»0.1.%» -5.5.0.1.6, .01 
2.23. 0.0.0. .5.0.0» 76.77 
124.1, 6.0.1.6. 6.0.1.1. .Ol 
49.2. 0.0.0, 0.0.-.5. 76.124 
225.1, 3.5.0..1» 3.5.0.0. .01 
224.1. 3.5.0..1. 3.75,0,.1, .01 


2.9. 0.0.0. .25,0.0. 225.226 

м 243,1, 6,0,.1, 6.0.0. .01 
21.1. 0.0.0, -9.5,0.0. 223.243 
265,1, -3.5.0,.1. -3.0..1, .01 
266.1. -3.0..1. -2.5.0..1» .01 
2.6. 0.0.0. 1.0.0. 265.26% 
278.001 


0.0.0.0. 20.557 





Figure 4.2 NEC Geometry Cards. 


A sample of the NEC data used to simulate the slot and 
the bulkhead is shown in Fig. 4.2. Sr 
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To produce the Numerical Green's Function (NGF) of this 
slot the data set shown in Appendix A was first runat a 
frequency of 20.357 MHz; the half-wavelength resonant 
frequency of the slot. 

It was mandatory to use the MVS (Multiple Virtual 
System) because the great majority of jobs were class C or G 
(CPU times per job of 3 minutes and 15 minutes , 
respectively). 


THETA = 60.00 ΡΗΙΞ 60.00 ETA = 90.00 





Figure 4.3 Wire Segment Used for Feeding the Slot. 


The next step was to feed the slot at different points 
along X-axis (from x=0.0Om through x=3.0m in 0.5m steps) апа 
to analyze its radiation patterns (vertical and horizontal) 
for E-field, its input impedance and to plot the near field 


inside the antenna. 
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Fig.4.3 shows the segment used to feed the structure 
oriented in the X-Y-Z coordinate system. 

Other data sets were developed to study the change of 
the slot input impedance with frequency variation. New 
NGF's had to be calculated for 5MHz, 10MHz and 30MHz. 

The slot width (w) was varied and the bandwidth (BW) 


performance studied. 


THETA = 60.00 PHI= 60.00 ЕТА = 90.00 





Figure 4.4 New Wire-Grid Model for w=0. 4m. 
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Figure 4.4 shows the new wire-grid model with w=0. 4m. 

After feeding the two structures at the same seven 
points as for w=0.2m, a comparison gave more information 
about the BW performance of these slot structures as "w" is 
varied. 

Values of input impedance were also included for the 
case where w-O. 5m. 

Keeping the same feed point (at x=0.0m), a comparison 
between the input impedance for "w" =0.2 and 0.4m, and four 
frequencies (5,10,20 and 30 MHz) was made. 

Finally, impedance is presented for frequencies between 
18 and 22 MHz, in 1 MHz steps, for both slots with w=0. 2m 
and w=0. 4m. 


30 


V. NUMERICAL RESULTS 


A.  INPUT IMPEDANCES 


TABLE I 
SLOT INPUT IMPEDANCE: F-20 MHZ AND W=0. 2M 


"x" (meters) Impedance (ohms) 


Distance 
464 - j 126.5 

443 = 5 113.7 
505 Уве >. 38 

300. | > 57) 

202. |. 0059 

TOT mor. 1 

34. mo: 





Seven data sets were run to evaluate the variation of 
input impedance of the slot as a function of "x", the 
feedpoint location along the horizontal axis. 

These results are shown on two different curves, one for 
resistance and the other for reactance, in Figures 5.1 and 
5.2. Table I contains the plotted values. 

Ihe point along x-axis where one would expect an input 
impedance of 50 ohms matched the theoretical value of 
Chapter II. In Fig.5.1, for example, it can be seen that for 
an input resistance of 50 ohms the value of "x" is close to 


2.88m, which matches the theory. One should notice that: 


К = 3.5 2 (A720) 2 3.5 = (300/20. 357/20) = 2.763m 
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Figure >. Input Resistance-20 Mhz. 


In Fig.5.2 it is easily seen that resonance is at a feed 
position of x-1. 8m. 

The same procedure was repeated for f=5 MHz as сап be 
seen in Appendix B. 

For 10 MHz and 30 MHz, the results appear respectively 
in Appendices C and D. 

Table II contains impedance results for w=0.2m and 


feeding distance of x=0. Om. 
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SLOT INPUT REACTANCE 
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Figure 5.2 Input Reactance-20 Mhz. 


Table III contains results for w=0.4m and feeding 
distance of x=0.0m. Here, the slot width is twice that of 


the previous case. 


ES 


TABLE te 
SLOT IMPEDANCE(X=0. OM AND ΝΞΟ. 2M):5 TO 30MHZ 


Frequency (MHz) Impedance( ohms) 
5549 0.028 + j 559 

1020 1.25 + j 142239 

2070 464.0 - 1 12635 

SONO 74.3 =- j 37% 


TABLE III 
SLOT IMPEDANCE(X=0. OM AND W=0. 4M):5 TO 30MHZ 


Frequency (MHz) Impedance( ohms) 
5. 0 0.253 E-3 +4 j 1115812 
ο 0.263 E-1 + 3 26”. 
20.0 πιο αυ - 1 792554 
3020 Бо - 3j 26739 





B. BANDWIDTH PERFORMANCE 

As explained in Chapter IV, the first task in this study 
was a comparison between both models (w=0.2m and w=0. 4m), 
keeping the frequency (f=20.0 MHz) constant. 

The results are in Fig. 5.3 and 5.4 where we see that 
for w=0.2m the input resistance and reactance both have 
smaller variation than those for  w-O. 4m. Thus there is an 
advantage in using the antenna with w=0.4m over that with 
W=0. 2m : a larger range of input impedances is available to 


WoOorR wich. 
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Figure 5.3 Slot Resistance Variation-20 MHz. 


Also included are two plots for w=0.5m, f=20MHz, one for 
the input resistance and the other for the input reactance, 
in Appendix E. 

Another possible parameter variation is with frequency, 
keeping the feed point location (x=0.0m) fixed, and studing 
Кие behavior of input impedance for the same slot width. 
This was done for 18 MHz to 22 MHz in steps of 1 MHz, for 
both cases where w-0.2m,x-0.0m (Table IV) and w=0. 4m, x=0. Om 
(Table V). 

Both plots corresponding to these tables are shown in 
ECCE El and Fig. Ε.2 It is easy to see that for the first 
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Figure 5.4 Slot Reactance Variation-20 MHz. 


value of w=0.2m the resonance is close to 19 MHz, while for 
w=0. 4m the resonance can be found ас #=20. 4 МЕ?! 
approximately. Thus one can say that the resonant frequency 


does not depend greatly on the slot width. 


С. NEAR FIELD INSIDE lee ore. 

The results of the near-field inside the slot for z=0. Om 
are shown in Fig. 5.5 

Referring to the geometry of Fig. 4.1, for 2=0.05т 
E-peak field is calculated just over the middle of the slot. 
The calculations меге repeated for z=0.0m and z=0.075m as 


shown in Appendix G. 
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TABLE IV 
SLOT IMPEDANCE(X=0. OM AND W=0. 2M): 18 TO 22MHZ 


Frequency (MHz) Impedance( ohms) 
1320 434.1 + j 279. 


19. 539.4 + j 125. 
20. 464.0 - j 126. 
one 373,3 Ро 
228 248.1 j 182. 


TABLE V 
SLOT IMPEDANCE(X=0. OM AND W=0. 4M): 18 TO 22MHZ 


Frequency (MHz) Impedance( ohms) 
18.0 22 мыл бы дену, 
τ; 56792 + j 4124.0 
20: 2105. 1 =$) 7928. / 
21 73618 Сат о 
22. 95.44 j 2546.0 





It was taken into account, for sake of simplification, 
that the field is symmetric with respect to the y-axis. 

Ihe E-peak is a composite of Ex and Ez components, with 
Ez being the dominant component. The smooth decline of 
E-peak along the x-direction can be observed in Fig.5.5 and 
Appendix G. 

There are some points of maximum along the near E-field 
curve as seen in Fig. 5.5. These points are exactly those 


where the junctions of the segments occur along the upper 
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REPRODUCED AT GOVERNMENT EXPENSE 
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Ејацке 5.5 Near E-Field 


and lower wires of the larger sides 
flat surface of the real bulkhead 


Wire grid. 


D. RADIATION PATTERNS 
1. Single Voltage Source 
The radiation pattern plots 
theory. 
For w=0.2m and £=20.0 MHz, 





for z=0. 050m. 


of the slot, because the 
is modeled here by the 


agreed with slot antenna 


two horizontal patterns 


(one for x=0.5m and the other for x=3,0m) and two vertical 


patterns (for the same feeding points along x-axis) were 


plotted. Appendix H contains these plots. 
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2. Multiple Voltage Sources 

One additional test with the same model (w=0.2m and 
f-20.0MHz) was considered. 

The slot pattern was not very sensitive to 
increasing the number of multiple feed segments. Appendix I 
shows this (horizontal patterns only) for three, nine and 
thirteen segments fed simultaneously. Very little change in 
shape is noted. 

Appendix J shows the NEC geometry card data set for 


multiple feeding using thirteen voltage sources. 
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VI. CONCLUSIONS AND SUGGESTIONS 


It has been shown that two main characteristics of a 
slot antenna were produced in this numerical modeling study. 
These were the antenna's terminal impedance and its 
radiation pattern. For the simulated slot antenna, the 
terminal impedance is the impedance seen by the voltage 
Source used to excite the antenna. This investigation 
produced good results and, as said before, this impedance is 
important with regard to matching the antenna to a 
transmission line and a source or receiver to assure 
efficiency in the transfer of power. The point along x-axis 
where one would expect an input impedance of 50 ohms matched 
the theoretical value of Chapter II. For "x" approximately 
equal to 2.88m, the input impedance was about 50 ohms. The 
resonance is at feed position of x=1.8m and, for that point, 
the input resistance is about 240 ohns. 

The feed point is not the only parameter affecting the 
computer model's impedance. The slot width variation led to 
a conclusion that there is an advantage in using antennas 
with broader "w" on account of having a larger range of 


input impedances to work with. 


Another possible parameter variation, keeping the feed 
point location fixed, is in frequency. Testing over the 
frequency range: 18 to 22 Miz; it was verified that the 


resonance does not depend greatly on the slot width. For 
both models with w=0.2m and w=0.4m, the resonant frequency 
was around 19 MHz and 20.4 MHz, respectively. 

The radiation pattern defines the directional gain of an 
antenna. In this thesis the spatial distribution of the phi- 
and theta-polarized electric field components radiated in 
the far-zone of antenna as well as near-fields were 


investigated. Both calculations agreed with the slot 
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antenna theory keeping іп mind that the group of points of 
maximum values appeared because а wire grid is only an 
approximation to a planar surface for currents on that 
surface. 

One suggestion is to place some type of dielectric 
backing in the slot region (Ref. 13], both for structural 
support and to seal against the environment. Therefore, it 
would be useful to develop a method for the wire-grid model 
to account for the presence of a slot dielectric. 

Another suggestion is to try some computation with 
boxed-in slot antennas and study their expected 
unidirectional radiation patterns. Another investigation is 
to perform a more detailed study of the influences of the 
slot width variation on bandwidth. 

In the same above reference, on page 7, Funke took into 
account, for sake of more accuracy,the "wire mesh density". 
This term actually deals with the length of the wire 
segments which form the grid for the modeled surface. He 
asked how tight must the mesh be in order to obtain 
reasonable results from the model. It must not be forgotten 
that the penalty for a tight mesh with very short wire 
segments is large computer storage area and long execution 
time. 

In this study, 556 segments were used (not including 
one required segment to represent the voltage source which 
feeds the slot antenna). Funke also began his investigation 
with a 555 segment model and showed that very little change 
in the calculated impedance was observed for 641 and 683 
segment wire-grid models. Decreasing the wire mesh density, 
the results were almost the same with a 383 segment model 
which would be a good start for further studies in this 
subject. 

From the results of Chapter IV and V it can be said that 


МЕС is an effective analysis tool for some slot antenna 


4] 


calculations, and should be further used in the application 


of low profile HF shipboard antenna designs. 
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APPENDIX A 


AN EXAMPLE OF JOB STREAM WITH NGF DATA SET 


//83 
//NEC 
//*NEC 

К GO 
//STEPLIB 
//FTO1FOO1 
//FTO4FOO1 
4 
//FTO5FOO1 
//FTO6FOO1 
//FTOSFOO1 
//ЕТ11Е001 


// 
//FT12F001 


Ж 
//ЕТ13Е001 


КА 
//FT14F001 


А 
//ЕТ15Е001 


// 
//ЕТ16Е001 


// 
//ЕТ20Е001 
2 
//ЕТ21Е001 


2 
//STEPNAME 


DAT 


JOB (1714,1234), 'NEC RUN',CLASS=G 


PROC VERSION=NPSNEC, STORAGE=1024K 
PROC VERSION=NPS1000, STORAGE=2048K 


EXEC PGM=&VERSION, REG I ON=&STORAGE 


DD 


DD 


DD 


DD 


DD 


DD 


DD 


DD 


DD 


DD 


DD 


DD 


DD 


DD 


DISP=SHR, DSN=MSS. F1595. NEC. LOAD 

DUMMY 

UNIT=SYSDA, SPACE=(CYL,(8,2)), 
DCB-( RECFM-VS , BLKSIZE-19069) 

DDNAME-SYSIN 

SYSOUT-* 

DUMMY 

UNIT=SYSDA, SPACE=(CYL,(8,2)), 
DCB-( RECFM-VS , BLKSIZE-19069) 

UNIT-SYSDA, SPACE-(CYL, (8,2) ), 
DCB=( RECFM=VS, BLKSIZE=19069) 

UNIT=SYSDA, SPACE=( CYL,(8,2)), 
DCB-( RECEM=VS, BLKSIZE=19069) 

UNIT=SYSDA, SPACE=( CYL,(8,2)), 
DCB=( RECEM=VS, BLKSIZE=19069) 

UNIT=SYSDA, SPACE=( CYL,(8,2)), 
DCB=( RECEM=VS, BLKSIZE=19069) 

UNIT=SYSDA, SPACE=(CYL,(8,2)), 
DCB=( RECFM=VS, BLKSIZE=19069) 

UNIT=SYSDA, SPACE-(CYL, (8,2) ), 
DCB-( RECFM=VS, BLKSIZE-19069) 

DUMMY 


PEND 
EXEC NEC,VERSION-DNPS1000,STORAGE-2048K 


Z/*DIRECT THE FILE 08 OUTPUT РАТА 
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ZZ/* TPO ITS КЕБЕК ЕЕЕ 
//*EITHER THE 5ТАСІМС DISK ОК М55 


Е 

ы 

//*GO. FTOBFOO1 DD DSN=S1714. S3. PLOTDATA, 

//* DISP=(NEW,KEEP), 

//* UNIT=3350, VOL=SER=MVSOO4, 

//* SPACE-(CYL,(4,4)), 

//* DCB=( LRECL=80, RECFM=FB, BLKSIZE=19040) 
ДЫ 

//*GO. FTOBFOO1 DD DSNZMSS. $1714. PLOTDATA( S3), 
А UNIT=3330V,VOL=SER=MS0402, 

//* SPAR ο. 

Ж DCB=( RECFM=FB, LRECL=80, 

a BLKSIZE=12960,DSORG=PO), 

//* DISP=( NEW, CATLG) 

Vine 

//*GO. ETO8FOO1 DD DSN=MSS. $1714. PLOTDATA(S3), 
ЖА DISP=(OLD,KEEP), 

//* DCB=( RECFM=FB, LRECL=80, BLKSIZE=12960, 
//* DSORG-PO) 

ДЕ 

//*GO. FTO8FOO1 DD DSN=S1714. S3. CURRDATA, 

Z: DISP=(NEW,KEEP), 

//* UNIT=3350,VOL=SER=MVS004, 

KR SPACE=(CYL,(4,4,), 

//* DCB=( LRECL=80,RECFM=FB, BLKSIZE=19040) 
Az 

//*GO. FTOBFOO1 DD DSN-MSS. $1714. CURRDATA( S3), 
АА UNIT=3330V,VOL=SER=MS0402, 

АЖ SPACE=( CYL,(4,4,2)), 

//* DCB=( RECEM=FB, LRECL=80,BLKSIZE=12960, 
и“ DSORG-PO), 

//* DISP-(NEW,CATLG) 

VE 
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//*GO. FTOBFOO1 DD DSN-MSS. 51714. CURRDATA( 53), 


2” DISP=(OLD,KEEP), 
Аз DCB=( RECFM=FB, LRECL=80, BLKSIZE=12960, 
5 DSORG=PO) 

е 


//*NEXT THE NGF FILE LOCATION. 
//*EITHER STAGING AREA DISK 
//*OR MSS DATA SET. 


у“ 

//*GO. FT20F001 DD DSN=MSS. $1714. NGFDATA( $3), 
//* DISP-SHR, 

//* DCB=( RECFM=VS, BLKSIZE=13030,DSORG=PO) 
Аш 

//*GO. FT20F001 DD DSN=MSS. $1714. NGFDATA( $3), 
//* DISP=(OLD,KEEP), 

//* DCB=( RECFM=VS, BLKSIZE=13030, DSORG=PO) 
Ы 

//%60. ЕТ20Е001 DD DSN-S1714. S3. NGFDATA, 

//* UNIT=3350, VOL=SER=MVSOO4, 

//* DISP=SHR 

Ж” 

//*G0. FT20F001 DD DSN=51714. FILE. FT20FOOL, 

//* UNIT=3350,VOL=SER=MVS004, 

//* DISP=SHR 

К" 


//* C CREATING THE NGF: 

//* C CRIANDO A MATRIZ COM A FUNCAO DE GREEN: 
//* C CRIANDO NOVA NGFDATA1 PARA CONTER S3 
//GO. FT20F001 DD DSN-MSS. 51714. ΝΟΕΡΑΤΑΙ( 53), 


иу UNIT=3330V,VOL=SER=MS0402, 

//* SPACE-(CYL,(4,4,2)), 

p DCB=( RECEM=VS, BLKSIZE=13030,DSORG=PO), 
DÅ DISP=( NEW, CATLG) 

ы 


//*GO. FT20FO01 DD DSN=51714. $3. NGFDATA, 


45 


yr DISP=( NEW, KEEP), 


Z UNIT=3350, VOL=SER=MVS004, 

а SPACE-(CYL,(4,4)), 

//* DCB=( RECFM-VS , BLKSIZE=19040) 
ы 

//*GO. ET20F001 DD DSN=S1714. FILE. FT20F001, 
//* DISP=(NEW, KEEP), 

Vie UNIT=3350 , VOL=SER=MVSOO4, 

E SPACE-(CYL,(4,4)), 

Vie DCB=( RECFM=VS, BLKSIZE=19040) 
ы 


//*NOW THE SOMMERFELD GRID LOCATION 

//% | 

//*GO. FT21FO01 DD DSNZ2MSS. $1714. SOMDATA( 53), 
//* DISP=SHR, 


ke: DCB=( RECFM=VS, BLKSIZE=13030, DSORG=PO) 
К 

//*GO. FT21FO001 рр DSN=S1714. S3. SOMDATA, 

2 UNIT=3350,VOL=SER=MVS004, 

τά. DISP-SHR 

i 


//*GO. FT21F001 DD DSN-8S1714. FILE. FT21FO001, 
UNIT=3350,VOL=SER=MVS004, 

//* DISP-SHR 

RE 

//GO. SYSIN DD * 

//* C INPUT DATA CARDS GO HERE 

CE BULKHEAD WITH SLOT ANTENNA 

GW 1,1, -6,0,4.6, -6,0,3.6, .01 

GW 2,1, -6,0,4:6, -5,0,496 (1 

cm 2,11 ο πο 

см 25,1, 6,0,4.6, 6,0,3.6, .01 

СМ: 25.20 060-0004 

GW 76,1, -6,0,1.6, -6,0,1.1, .01 
GW 77,1, -6,0,1.6, -5.5,0,1.6, .01 
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GM 
GW 
GM 
GW 
GW 
GM 
GW 
GM 
GW 
GW 
GM 
GX 
GE 
FR 
WG 
EN 


Hie 


/* 


το О О О, 0,0, 76.77 
meme 016, 6,0 1.1,.,.01 
H2 10 09.0.:0,.0,2957 75. 124 
ο ο GE ο ο ο, 791 
22 n 5011137750, . 1, .01 
Zu ШОО, Ол 25,0,050223. 224 
2495115. 6,07. 17.755040, 4 ВН 
Ziel. ο ο ο. 99. 5050702237243 
ZONE. -3.9,0,.l1. -Ә 6.1, 201 
2.5, τ, ον 03 ορ ο ο αἱ 
2ο ο ο ο -αἶ 00572657 266 
278,001 


ο. 96, 5.0 
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APPENDIX B 
SLOT INPUT IMPEDANCE FOR F-5MHZ AND W=0. 2M 


u HH 





Figure B.1 Input Resistance-f-5MHz and w-O. 2m. 
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Figure B.2 Input Reactance-f=5MHz and w=0.2m. 
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APPENDIX C 
SLOT INPUT IMPEDANCE FOR F=10MHZ AND W=0. 2M 
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Figure C.1 Input Resistance-f-10MHz and w=0. Zm. 
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Figure C.2 Input Reactance-f-10MHz and w=0. 2m. 
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REPRODUCED AT GOVERNMENT EXPENSE 


APPENDIX D 
SLOT INPUT IMPEDANCE FOR F-30MHZ AND W-O. 2M 


SLOT INPUT RESISTANCE 


211111. а 
1111 1 1 ee 
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сте рт Input Resistance-f=30MHZ and w=0. Zm. 
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Figure D. 2 Input Reactance-f=30MHz and w=0. 2m. 


REPRODUCED AT GOVERNMENT EXPENSE 
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APPENDIX E 
SLOT INPUT IMPEDANCE FOR F=20MHZ AND W=0. 5M 


SLOT INPUT RESIS TEAOR 
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Figure E. 1 Input Resistance-f-20MHz and w-O. 5m. 
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Slot Input Impedance-w=0. 2m and x=0. Om. 
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SLOT INPUT IMPEDANCE 
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Figure Е.2 Slot Input Impedance-w=0. 4m and x=0. Om. 
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APPENDIX 6 
NEAR FIELD INSIDE ΤΗΕ Τη. 
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Figure G.1 Near-Field for z=0. Om. 
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NEAR E-FIELD INSIDE SLOT 
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Figure G.2 Near-Field for z=0. 075m. 
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SEGMENTS FED:THIRTEEN 





Figure I.3 E-Field Horizontal Pattern : 13 sources. 
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PPENDIX J 


DATA SET FOR A MULTIPLY FED SLOT ANTENNA 


//GO. SYSIN DD * 

CM FEEDING THE SLOT AT THIRTEEN POINTS 
CE GETTING THE NGF FOR w=0.2m and £=20.0MHz 
СЕ 

ο σος τ ο ο. 100. 0. , 0. ;-. 100,. 01 
см 601,1, .5,0,.100, .5,0,-.100,.01 
GW 602,1, -.5,0,.100, -.5,0,-.100,.01 
GW 603,1, -1.,0,.100, -1.,0,-.100,.01 
ο ει ІП о, ТОО, 1. ,0,-. 100,. 01 
GW 605,1, 1.5,0,.100, 1.5,0,-.100,.01 
GW 606,1,-1.5,0,.100,-1.5,0,-.100,.01 
GW 607,1,-2.0,0,.100,-2.0,0,-.100,.01 
en ϱο 250, 1000 2.0,0,-.100,.01 
οὐ. 25,0, 100, 2.5,0,-. 100,. 01 
GW 610,1,-2.5,0,.100,-2.5,0,-.100,.01 
GN on 7-3 о. 0. 100.-3.0.0,-.100,.01 
ЕСЕ 00 19907 3.0,0,-.109,.01 
СЕ 

ЕХ 0,600,1 

ЕХ 0,601,1 

ЕХ 0,602,1 

ЕХ 0,603,1 

ЕХ 0,604,1 

ЕХ 0,605,1 

ЕХ 0,606,1 

ЕХ 0,607,1 

EX 0,608,1 

EX 0,609,1 

EX 0,610,1 

EX 0,611,1 
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